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Using SU(3) lattice QCD, we perform the detailed studies of the three-quark 
and the multi-quark potentials. From the accurate calculation for more than 
300 different patterns of 3Q systems, the static ground-state 3Q potential Kjq S ' 
is found to be well described by the Coulomb plus Y-type linear potential (Y- 
Ansatz) within 1%-level deviation. As a clear evidence for Y-Ansatz, Y-type 
flux-tube formation is actually observed on the lattice in maximally- Abelian 
projected QCD. For about 100 patterns of 3Q systems, we perform the accurate 
calculation for the 1st excited-state 3Q potential V3Q 8 ' by diagonalizing the QCD 
Hamiltonian in the presence of three quarks, and find a large gluonic-excitation 
energy A_E 3 q = V£q — V$q' of about 1 GeV, which gives a physical reason 
of the success of the quark model. AE 3 q is found to be reproduced by the 
"inverse Mercedes Ansatz", which indicates a complicated bulk excitation for 
the gluonic-excitation mode. We study also the tetra-quark and the penta- 
quark potentials in lattice QCD, and find that they are well described by the 
OGE Coulomb plus multi-Y type linear potential, which supports the flux-tube 
picture even for the multi-quarks. 



1 Introduction 

In 1966, Yoichiro Nambu [I] first proposed the SU(3) gauge theory, i.e., quantum 
chromodynamics (QCD), as a candidate for the fundamental theory of the strong 
interaction, just after the introduction of the "new" quantum number, "color" 

2 . In 1973, the asymptotic freedom of QCD was theoretically proven [3], and, 
through the applicability check of perturbative QCD to high-energy hadron 

1 Invited Talk presented by Hideo Suganuma at the International Conference on "High 
Energy and Mathematical Physics" (ICHEMP05), 4-7 April 2005, Cadi Ayyad University, 
Marrakech, Morocco. 
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reactions, QCD has been established as the fundamental theory of the strong 
interaction. Even at present, however, it is very difficult to deal with QCD 
analytically due to its strong-coupling nature in the infrared region. Indeed, in 
spite of its simple form, QCD creates thousands of hadrons and leads to various 
interesting nonperturbative phenomena such as color confinement ^ E] and 
dynamical chiral-symmetry breaking |7j. Instead, lattice QCD has been applied 
as the direct numerical analysis for nonperturbative QCD. 

In 1979, the hrst application [5] of lattice QCD Monte Carlo simulations was 
done for the inter-quark potential between a quark and an antiquark using the 
Wilson loop. Since then, the study of the inter-quark force has been one of the 
important issues in lattice QCD [§]. Actually, in hadron physics, the inter-quark 
force can be regarded as an elementary quantity to connect the "quark world" 
to the "hadron world", and plays an important role to hadron properties. 

In this paper, we perform the detailed and high-precision analyses for the 
inter-quark forces in the three-quark and the multi-quark systems with SU(3) 
lattice QCD [TOl EH H21 H31 C3I CHI CHI 113 HBl HSl III], and try to extract the 
proper picture of hadrons. 

2 The Three-Quark Potential in Lattice QCD 

In general, the three-body force is regarded as a residual interaction in most 
fields in physics. In QCD, however, the three-body force among three quarks is 
a "primary" force reflecting the SU(3) gauge symmetry. In fact, the three-quark 
(3Q) potential is directly responsible for the structure and properties of baryons, 
similar to the relevant role of the QQ potential for meson properties. Further- 
more, the 3Q potential is the key quantity to clarify the quark confinement in 
baryons. However, in contrast to the QQ potential [H], there were almost no 
lattice QCD studies for the 3Q potential before our study in 1999 P0]> in spite 
of its importance in hadron physics. 

As for the functional form of the inter-quark potential, we note two theoret- 
ical arguments at short and long distance limits. 

1. At the short distance, perturbative QCD is applicable, and therefore inter- 
quark potential is expressed as the sum of the two-body OGE Coulomb 
potential. 

2. At the long distance, the strong-coupling expansion of QCD is plausible, 
and it leads to the flux-tube picture |21j . 

Then, we theoretically conjecture the functional form of the inter-quark poten- 
tial as the sum of OGE Coulomb potentials and the linear potential based on the 
flux-tube picture. Of course, it is highly nontrivial that these simple arguments 
on UV and IR limits of QCD hold for the intermediate region. Nevertheless, 
the QQ potential Vqq(r) is well described with this form as [Ql llllIT^] 

V QQ (r) = -^+<J QQ r + C QQ . (1) 
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For the 3Q system, there appears a junction which connects the three flux- 
tubes from the three quarks, and Y-type flux-tube formation is expected ^] 
E3I2JI221I2B1- Therefore, the (ground-state) 3Q potential is expected to be the 
Coulomb plus Y-type linear potential, i.e., Y-Ansatz, 

^3Q ' = -^3Q £ T-^-T + <J3QL min + C 3Q , (2) 

where L m [ n is the length of the Y-shaped flux-tube. 

For more than 300 different patterns of spatially-fixed 3Q systems, we cal- 
culate the ground-state 3Q potential V^q' from the 3Q Wilson loop W3Q using 
SU(3) lattice QCD [TT1 H21 HS1 HI with the standard plaquette action at the 
quenched level on various lattices, i.e., (/3=5.7, f2 3 x 24), (/3=5.8, 16 3 x 32), 
(,0=6.0, 16 3 x 32) and ((3 = 6.2, 24 4 ). For the accurate measurement, we con- 
struct the ground-state-dominant 3Q operator using the smearing method. Note 
that the lattice QCD calculation is completely independent of any Ansatz for 
the potential form. 

To conclude, we find that the static ground-state 3Q potential V^q S ' is well 
described by the Coulomb plus Y-type linear potential (Y-Ansatz) within 1%- 
level deviation |11M12| . To demonstrate this, we show in Fig.l the 3Q confine- 
ment potential V^Q ni , i.e., the 3Q potential subtracted by the Coulomb part, 
plotted against the Y-shaped flux-tube length L mm . For each (3, clear linear 
correspondence is found between the 3Q confinement potential VJjq and L m [ n , 
which indicates Y-Ansatz for the 3Q potential. 
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Figure 1: The 3Q confinement potential V 3 Q n , i.e., the 3Q potential subtracted 
by the Coulomb part, plotted against the total flux-tube length L m ; n of Y-Ansatz 
in the lattice unit. 
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Figure 2: The lattice QCD result for Y-type flux-tube formation in the spatially- 
hxed 3Q system in MA projected QCD. The distance between the junction and 
each quark is about 0.5 fm. 



Recently, as a clear evidence for Y-Ansatz, Y-type flux-tube formation is 
actually observed in maximally-Abelian (MA) projected lattice QCD from the 
measurement of the action density in the spatially-fixed 3Q system ^3 ^3 ^3 
I24| . (See Fig. 2.) In this way, together with recent several analytical studies 
and other lattice-QCD work [23, Y-Ansatz for the static 3Q potential 
seems to be almost settled. 



3 The Gluonic Excitation in the 3Q System 

In this section, we study the excited-state 3Q potential and the gluonic ex- 
citation in the 3Q system using lattice QCD E|. The excited-state 3Q 
potential V 3 e q ' is the energy of the excited state in the static 3Q system. The 
energy difference AE^q = V£q — V^q' between V^q' and V£q is called as the 
gluonic-excitation energy, and physically means the excitation energy of the 
gluon-field configuration in the static 3Q system. In hadron physics, the glu- 
onic excitation ^UEIEEI is one of the interesting phenomena beyond the quark 
model, and relates to the hybrid hadrons such as qqG and qqqG in the valence 
picture. 

For about 100 different patterns of 3Q systems, we calculate the excited-state 
potential in SU(3) lattice QCD with 16 3 x 32 at (3=5.8 and 6.0 at the quenched 
level by diagonalizing the QCD Hamiltonian in the presence of three quarks. 
In Fig. 3, we show the 1st excited-state 3Q potential V 3 q s ■ and the ground-state 
potential Vjq 3 '. The gluonic excitation energy AE^q = V£q' — V$q' in the 3Q 
system is found to be about lGeV in the hadronic scale as 0.5fm < L min < 
1.5fm. Note that the gluonic excitation energy of about lGeV is rather large 
compared with the excitation energies of the quark origin. This result predicts 
that the lowest hybrid baryon qqqG has a large mass of about 2 GeV. 
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Figure 3: The 1st excited-state 3Q potential V£q and the ground-state 3Q 
potential V 3 q S '. The lattice results at (3 — 5.8 and j3 — 6.0 well coincide except 



for an irrelevant overall constant. The gluonic excitation energy A.E3Q 
V 3 q S ' is about lGeV in the hadronic scale as 0.51m < L m i n < 1.51m. 



T/C.S. 



3.1 Inverse Mercedes Ansatz for the Gluonic Excitation 

Next, we investigate the lunctional form of the gluonic excitation energy in 3Q 
systems AE 3 q = V 3 c q- — V^q' , where the Coulomb part is expected to 

be canceled between V 3 q' and Kjg'. After some trials, as shown in Fig. 4, we 
find that the lattice data of the gluonic excitation energy AE 3 q = V 3 c q' — V 3 6 q' 
are relatively well reproduced by the "inverse Mercedes Ansatz" [T2|, AE 3 q = 

^ + g, Ly = y,U V*? - + e 2 = \ E i#i p^Q7 ^pq 4 , ^pp,), 

where L Y denotes the "modified Y-length" defined by the half perimeter of the 
"Mercedes form" as shown in Fig.4(a). As for (K, G, £), we find (K ~ 1.43, 
G ~ 0.77 GeV, £ - 0.116 fm) at /3 = 5.8, and [K ~ 1.35, G ~ 0.85 GeV, £ ~ 
0.103 fm) at /3 = 6.0. 




Figure 4: (a) The Mercedes form for the 3Q system, (b) & (c) Lattice QCD 
results of the gluonic excitation energy A^q = V^q 8 ' — V 3 S q' in the 3Q system 
plotted against the modified Y-length Ly. The dashed curve denotes the inverse 
Mercedes Ansatz. 



The inverse Mercedes Ansatz indicates that the gluonic-excitation mode is 
realized as a complicated bulk excitation of the whole 3Q system. 
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3.2 Behind the Success of the Quark Model 



Here, we consider the connection between QCD and the quark model in terms of 
the gluonic excitation ^3EJE3E1- While QCD is described with quarks and 
gluons, the simple quark model successfully describes low-lying hadrons even 
without explicit gluonic modes. In fact, the gluonic excitation seems invisible 
in low-lying hadron spectra, which is rather mysterious. 

On this point, we find the gluonic-excitation energy to be about lGeV or 
more, which is rather large compared with the excitation energies of the quark 
origin. Therefore, the contribution of gluonic excitations is considered to be neg- 
ligible and the dominant contribution is brought by quark dynamics such as the 
spin-orbit interaction for low-lying hadrons. Thus, the large gluonic-excitation 
energy of about lGeV gives the physical reason for the invisible gluonic excita- 
tion in low-lying hadrons, which would play the key role for the success of the 
quark model without gluonic modes ^3 G3 EH • 

In Fig. 5, we present a possible scenario from QCD to the massive quark 
model in terms of color confinement and dynamical chiral-symmetry breaking. 

Quantum Chromodynamics 



(Strong) Color Confinement 



(Tightly Stretched) 
Color Flux- Tube Formation 
with a Large String Tension 
a ~ 1 GeV/fm 



Large Excitation Energy 
of the Flux- Tube Vibration 



Large Gluonic Excitation 
Energy AE ~ 1 GeV 



Absence of Gluonic Excitation 
Modes in Low-lying Hadrons 



Dynamical Chiral Symmetry 
Breaking 



Larf 



i Constituent Quark Mass 
M, ~ 300 MeV 



Non-Relativistic Treatment 
on Quark Dynamics 
in Low-lying Hadrons 



Massive Quark Model for Low-lying Hadrons 



Figure 5: A possible scenario from QCD to the quark model in terms of color 
confinement and dynamical chiral-symmetry breaking (DCSB). DCSB leads to 
a large constituent quark mass of about 300 MeV, which enables the non- 
relativistic treatment for quark dynamics approximately. Color confinement 
results in the color flux-tube formation among quarks with a large string tension 
of a ~ 1 GcV/fm. In the flux-tube picture, the gluonic excitation is described 
as the flux-tube vibration, and its energy is expected to be large in the hadronic 
scale. The large gluonic-excitation energy of about 1 GeV leads to the absence 
of the gluonic mode in low-lying hadrons, which plays the key role to the success 
of the quark model without gluonic excitation modes. 
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4 Lattice QCD Study for the Inter-Quark Inter- 
action in Multi-Quark Systems 

In this section, we perform the first study of the inter-quark interaction in multi- 
quark systems in lattice QCD [H [T7| El HH1 1201 , and seek for the QCD-based 
quark-model Hamiltonian to describe multi-quark hadrons |29) . The quark- 
model Hamiltonian consists of the kinetic term and the potential term, which 
is not known form QCD in multi-quark systems. 

As for the potential at short distances, the perturbative one-gluon-exchange 
(OGE) potential would be appropriate, due to the asymptotic nature of QCD. 
For the long-range part, however, there appears the confinement potential as a 
typical non-perturbative property of QCD, and its form is highly nontrivial in 
the multi-quark system. 

In fact, to clarify the confinement force in multi-quark systems is one of the 
essential points for the construction of the QCD-based quark-model Hamilto- 
nian. Then, in this paper, we investigate the multi-quark potential in lattice 
QCD, with paying attention to the confinement force in multi-quark hadrons. 

4.1 The OGE Coulomb plus Multi-Y Ansatz 

We first consider the theoretical form of the multi-quark potential, since we will 
have to analyze the lattice QCD data by comparing them with some theoretical 
Ansatz. 

By generalizing the lattice QCD result of the Y- Ansatz for the three-quark 
potential, we propose the one-gluon-exchange (OGE) Coulomb plus multi-Y 
Ansatz [Hil fT71 1T%1 fT31 1277) . 

2 rparpa 

for the potential form of the multi-quark system. Here, the confinement po- 
tential is proportional to the minimal total length L m ; n of the color flux tube 
linking the quarks, which is multi-Y shaped in most cases. 

In the following, we study the inter-quark interaction in multi-quark systems 
in lattice QCD, and compare the lattice QCD data with the theoretical form in 
Eq.©. Note here that the lattice QCD data arc meaningful as primary data 
on the multi-quark system directly based on QCD, and do not depend on any 
theoretical Ansatz. 

4.2 Formalism of the Multi-Quark Wilson Loop 

Next, we formulate the multi-quark Wilson loop to obtain the multi-quark po- 
tential in lattice QCD EH [T7| GJ 033 EDI . 

Similar to the derivation of the QQ potential from the Wilson loop, the 
static multi-quark potential can be derived from the corresponding multi-quark 
Wilson loop. We construct the tetraquark Wilson loop Waq and the pentaquark 
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Wilson loop W$q in a gauge invariant manner as shown in Figs. 6(a) and (b), 
respectively. 




Figure 6: (a) The tetraquark Wilson loop W^q. (b) The pentaquark Wilson 
loop W$q. The contours M, Mi, Rj, Lj(i = 1,2, j = 3,4) are line-like and 
Rj,Lj(j = 1,2) are staple-like. The multi-quark Wilson loop physically means 
that a gauge-invariant multi-quark state is generated at t — and annihilated 
at t = T with quarks being spatially fixed in R 3 for < t < T. 

The tetraquark Wilson loop W4Q and the pentaquark Wilson loop W5Q are 
defined by 

W 4 q = ^tr(Mifii 2 M 2 ii2), 

W 5Q = ie afcc e a ' b ' c 'M aa '(i? 3 i? 12 i? 4 )'"''(i 3 ii2i4) cc ', (4) 
where M, Mi, Lj and Rj (i=l,2, j=l, 2,3,4) are given by 

M, Mi, Rj,Lj = Pcxp{ig [ dx^A^x)} G SU(3) C . (5) 

JM,Mi,Rj,Lj 

Here, R12 and Li 2 arc defined by 

p£ a ' a — }_^abc ^a'b' c' j^bb' j^cc' jjx 'a _ ^_^abc^a'b' c jbb' j^cc' /g\ 
12 ^ 12' 12 ^ 12' V / 

The multi-quark Wilson loop physically means that a gauge-invariant multi- 
quark state is generated at t — and annihilated at t = T with quarks being 
spatially fixed in R 3 for < t < T. 

The multi-quark potential is obtained from the vacuum expectation value of 
the multi-quark Wilson loop as 

V iQ = lim hn(W 4Q ), V 5Q = lim hn(W 5Q ). (7) 

1 —*oo ± 1 — *oo ± 

4.3 Lattice QCD Result of the Pentaquark Potential 

We perform the first study of the pentaquark potential V5Q in lattice QCD 
with (/3=6.0, 16 3 x 32) for 56 different patterns of QQ-Q-QQ type pentaquark 
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configurations, as shown in Fig. 7. As the conclusion, the lattice QCD data of 
^5Q are found to be well described by the OGE Coulomb plus multi-Y Ansatz, 
i.e., the sum of the OGE Coulomb term and the multi-Y-type linear term based 
on the flux-tube picture [15111711151 1131 120] . 



3 3 




Figure 7: A QQ-Q-QQ type pentaquark configuration. In the 5Q system, 
(QijQa) and (Q3,Q4) form 3 representation of SU(3) color, respectively. The 
lattice QCD results indicate the multi- Y-shaped flux-tube formation in the QQ- 
Q-QQ system. 



We show in Fig. 8 the lattice QCD results of the 5Q potential Ve>q for sym- 
metric planar 5Q configurations as shown in Fig. 7, where each 5Q system is 
labeled by d = Q1Q2/2 and h = Q1Q3. 
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Figure 8: Lattice QCD results of the pentaquark potential V5Q for symmetric 
planar 5Q configurations in the lattice unit: (a) V$q v.s. (d, h) and (b) V5Q v.s. 
Zvmkr The symbols denote the lattice QCD data, and the curves the theoretical 
form of the OGE plus multi-Y Ansatz. 

In Fig. 8, we add the theoretical curves of the OGE Coulomb plus multi-Y 
Ansatz, where the coefficients (j45Q,<75q) are set to be (^43Q,CT3q) obtained 
from the 3Q potential ^2J. (Note that there is no adjustable parameter in the 
theoretical Ansatz apart from an irrelevant constant.) In Fig. 8, one finds a good 
agreement between the lattice QCD data of V5Q and the theoretical curves of 
the OGE Coulomb plus multi-Y Ansatz. 

In this way, the pentaquark potential V5Q is found to be well described by 
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the OGE Coulomb plus multi-Y Ansatz as [T5I IT71IT31ITg| 

V 5Q = - A 5Q {(— + _) + _(_ + _ + _ + _ ) 

ri2 ?"34 2 n 5 r 25 r 35 r 45 

where r%j is the distance between Qi and Qj. Here, L^ n is the minimal total 
length of the flux tube, which is multi- Y-shaped as shown in Fig. 7. This lattice 
result supports the flux-tube picture for the 5Q system. 



4.4 Tetraquark Potential and Flip-Flop in Lattice QCD 

For about 200 different patterns of QQ-QQ configurations, we perform the de- 
tailed study of the tetraquark potential Vaq in lattice QCD with /3=6.0, 16 3 x 32, 
and find the following results [T?)l HSI . 

1. When QQ and QQ are well separated, the 4Q potential Vaq is well de- 
scribed by the OGE Coulomb plus multi-Y Ansatz, which indicates the 
multi- Y-shaped flux-tube formation as shown in Fig. 9 (a). 

2. When the nearest quark and antiquark pair is spatially close, the 4Q 
potential Vaq is well described by the sum of two QQ potentials, which 
indicates a "two-meson" state as shown in Fig. 9(b). 




Figure 9: (a) A connected tetraquark (QQ-QQ) configuration and (b) A "two- 
meson" configuration. The lattice QCD results indicate the multi- Y-shaped 
flux-tube formation for the connected 4Q system. 



We show in Fig. 10 the lattice QCD results of the 4Q potential V4Q for 
symmetric planar 4Q configurations as shown in Fig. 9, where each 4Q system 
is labeled by d = Q1Q2/2 and h = Q1Q3. 

For large value of h compared with d, the lattice data seem to coincide with 
the solid curve of the OGE Coulomb plus multi-Y Ansatz, 
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Figure 10: Lattice QCD results of the tetraquark potential V4Q for symmetric 
planar 4Q configurations in the lattice unit. The symbols denote the lattice 
QCD data. The solid curve denotes the OGE plus multi-Y Ansatz, and the 
dottcd-dashcd curve the two-meson Ansatz. 



where Lj^ is the minimal total length of the flux tube, which is multi- Y-shaped 
as shown in Fig. 9(a). Here, the coefficients (A4Q, <J4q) are set to be (^3Q, 03q) 
obtained from the 3Q potential |12| . 

For small h, the lattice data tend to agree with the dotted-dashed curve of 
the "two-meson" Ansatz, where the 4Q potential is described by the sum of two 
QQ potentials as V Q Q(r 13 ) + Vqq^) = 2V q q(}i). 

Thus, the tetraquark potential V4Q is found to take the smaller energy of 
the connected 4Q state or the two-meson state. In other words, we observe a 
clear lattice QCD evidence of the "flip-flop", i.e., the flux-tube recombination 
between the connected 4Q state and the two-meson state. This lattice result 
also supports the flux-tube picture for the 4Q system. 



4.5 Proper Quark-Model Hamiltonian for Multi-Quarks 

From a series of our lattice QCD studies [1111111131111111111111311111111101 

on the inter-quark potentials, the inter-quark potential is clarified to consist 
of the one-gluon-exchange (OGE) Coulomb part and the flux-tube-type linear 
confinement part in QQ-mesons, 3Q-baryons and multi-quark (4Q, 5Q) hadrons. 

Furthermore, from the comparison among the QQ, 3Q, 4Q and 5Q potentials 
in lattice QCD, we find the universality of the string tension a, 

a QQ — °3Q — O4Q ~ <T5Q, (10) 

and the OGE result of the Coulomb coefficient A as 

2 A QQ ~ A 3Q ~ A *Q ~ A *Q ( n ) 

in Eqs.Q, @, © and @. 

Here, the OGE Coulomb term is considered to originate from the OGE 
process, which plays the dominant role at short distances, where perturbative 
QCD is applicable. The flux-tube-type linear confinement would be physically 
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interpreted by the flux-tube picture, where quarks and antiquarks are linked by 
the one-dimensional squeezed color-electric flux tube with the string tension a. 

To conclude, the inter-quark interaction would be generally described by 
the sum of the short-distance two-body OGE part and the long-distance flux- 
tube- type linear confinement part with the universal string tension a ~ 0.89 
GeV/fm. 

Thus, based on the lattice QCD results, we propose the proper quark-model 
Hamiltonian H for multi-quark hadrons as 

^ = 1^+ M i +J2 V OGE+ ^min, (12) 

i i<j 

where L m i n is the minimal total length of the flux tube linking quarks. Vq GE 
denotes the OGE potential between ith and jth quarks, which becomes the 
OGE Coulomb potential in Eq.© in the static case. Mj denotes the constituent 
quark mass. The semi-relativistic treatment would be necessary for light quark 
systems. 

It is desired to investigate various properties of multi-quark hadrons with 
this QCD-based quark model Hamiltonian H . 
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